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Abstract 

Aims. In this paper we explore the possibility that the recently detected reflected light signal of 51 Peg b could be caused by a ring 
system around the planet. 

Methods. We use a simple model to compare the observed signal with the expected signal from a short-period giant planet with rings. 
We also use simple dynamical arguments to understand the possible geometry of such a system. 

Results. We provide evidence that, to a good approximation, the observations are compatible with the signal expected from a ringed 
planet, assuming that the rings are non-coplanar with the orbital plane. However, based on dynamical arguments, we also show that 
this configuration is unlikely. In the case of coplanar rings we then demonstrate that the incident flux on the ring surface is about 2% 
the value received by the planet, a value that renders the ring explanation unlikely. 

Conclusions. The results suggest that the signal observed cannot in principle be explained by a planet-l-ring system. We discuss, 
however, the possibility of using reflected light spectra to detect and characterize the presence of rings around short-period planets. 
Finally, we show that ring systems could have already been detected by photometric transit campaigns, but their signal could have 
been easily misinterpreted by the expected light curve of an eclipsing binary. 

Key words. (Stars:) Planetary systems. Planets and satellites: detection. Techniques: spectroscopy. Stars: Individual: 51 Peg 


1. Introduction 

The detection of the atmospheres of extrasolar planets is becom¬ 
ing one of the major research topics in the exoplanet field (for 
a recent review see Burrows 2014 1 . Current technology and a 


detailed data analysis have already allowed the signature of the 
atmospheres of other worlds to be detected using different meth¬ 
ods, such as transmission spectroscopy (e.g., Charbonneau et al. 


[2002||Vidal-Madjar et al.|200^ [Madhusudhan et al.|2014|l. 


oc- 


cultations (e.g., |Deming et al.||2005[ [Demory et al. [2012| l, and 

phase curve variations (e.g., Angerhausen et al. 2()14[ ). These 
studies allowed several detailed analyses of exoplanet atmo¬ 
spheres, including tracing of thermal or albedo maps of the plan¬ 
ets (e.g., |Knutson et al. 2007 1 |Stevenson et al.||2014[ |Demory 
let al.|2013| l. 

Although a large majority of the exoplanet atmosphere stud¬ 
ies involved space-based data, the use of ground-based instru¬ 


mentation to detect exoplanet atmospheres is providing a grow¬ 
ing amount of information. This is particularly true concerning 
the use of high-resolution spectroscopic techniques. Using both 
optical and the near-infrared (near-IR) wavelengths, these meth¬ 
ods allowed the spectrum to be probed in detail for several exo¬ 
planets (for some examples see|Snellen et al.|2010||Birkby et al. 
20T3j|Wyttenbach et al.|20T5]l. 


In a recent paper, Martins et al. ( 2015| ) have explored a new 
technique for detecting the signature of a high-resolution (op¬ 
tical) reflected light spectrum from an exoplanet. This detec¬ 
tion allowed estimation of the radius and albedo of the histori¬ 
cal 51 Pegb planet (Mayor & Queloz 1995|l, suggesting that this 


planet may be a high-albedo, inflated 
7 b (with Ag = 0.35 Demory et al.|20lJ 


lot-Jupiter such as Kepler- 


The predicted star-to-planet flux ratio for a star-i-planet sys¬ 


tem be estimated from (e.g., Seager|2010|l: 




( 1 ) 


* Based on observations collected at ESO facilities under program 
091.C-0271 (with the HARPS spectrograph at the ESO 3.6-m telescope. 
La Silla-Paranal Observatory). 


where Ag is the geometric albedo of the planet, a the semi¬ 
major axis of the orbit, g{a) the phase function, and the 
planetary radius. The FWHM (22.6+3.6 km s ') and amplitude 
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(6.0+0.4x10“^) of the detected planet-cross-correlation function 
(CCF), as detected in [Martins et al. (2015 i, when compared 
to the values of the stellar CCF (7.47kms“‘ and 0.48, respec¬ 
tively), would lead to a planet-to-star flux ratio of 3.8x10 By 
applying the equation above, we would then derive a geomet¬ 
ric albedo far above unity if we assume a jovian-like radius for 
Sl PeglQ 

[Martins et al. ( |2015| l presented the results of some simula¬ 
tions suggesting that the observed (and larger than expected) 
FWHM broadening can be an artifact produced by non-Gaussian 
noise in the data, together with the fact that their detection was 
only possible at a three-sigma level. The authors thus only used 
a comparison of the CCF depths to derive indicative values for 
the albedo and planetary radius. (The two parameters are degen¬ 
erate.) Martins et al. also suggested that the parameters of the 
CCF should be taken just as indicative, even if they consider the 
detection solid. 

It is, however, interesting to explore the possibility that the 
observed CCF values are real. In this case, what could explain 
such a wide and deep CCF as observed? One possibility for ex¬ 
plaining the large FWHM would be the presence of strong winds 
or a very fast rotation velocity (close to the observed FWHM). 
The signature of strong winds in exoplanets has indeed been 
observed using transmission spectroscopy (e.g., [Snellen et al.| 
2010[l. However, such broadening of the CCF would also im¬ 


ply a decrease in its depth. The wind explanation would thus 
not be able to explain the total surface of the CCF. Some extra 
component would be needed. Besides this, though not discussed 
directly in the paper, the near-IR signal of 51 Peg b detected by 
[Brogi et ar] ( |2013[ l using CO lines does not seem to show a clear 
sign of extra broadening (see their Figure 3). Winds should pro¬ 
duce a broadening that is independent of the wavelength of the 
observations. 

In the present paper we explore a new interpretation of the 
measurements done by Martins et al. In particular, we try to un¬ 
derstand if the signal detected could be explained by a reflective 
ring system around 51 Pegb. The existence of such rings has 
already been demonstrated to be possible around hot-Jupiters 
(e.g., Schlichting & Chang[[201 1| . In Sect. [^ we unv eil our hy¬ 
pothesis based on the observations of |Martins et al. (2015 i, and 
we explore how a ring system should be able to explain the ob¬ 
served signal at opposition, i.e., when the Earth, the star, and 
the planet are almost aligned (in the same order). In this study, 
we assume that the rings are not coplanar with the planet or¬ 
bit and provide the results as a function of the angle <p between 
the Earth-star-planet line and the ring plane. In Sect.we use 
dynamical constraints to investigate whether the necessary con¬ 
figuration corresponds to a physical scenario. A more detailed 
model of the reflected light from a coplanar ring system is then 
explored in Sect. [3.1[ In Sect.we briefly show that rings may 
have already been detected using transit photometry, though in 
some cases their signature could have been interpreted as caused 
by the eclipse of a stellar companion. We conclude in Sect.|^ 


2. The case for rings: a simple model 

In Eq.[T]we denote the planet-to-star flux ratio expected for a 
star-(-planet system. The presence of rings around the planet 
would alter this ratio, because these would also reflect light to¬ 
ward the observer. 


* Note that the area of a Gaussian is proportional to 
FWHMx Amplitude 


The orbital inclination of 51 Peg b is close to 80 degrees 
( Brogi etaL][2013 Martins et al. 2()15[ ). Moreover, in Martins 
et al.[ ( [2015[ ), the detection of the reflected light has been done 
almost at opposition. Thus, to get a rough estimation of the re¬ 
flected light in this condition, we assume a simple geometry of 
the problem in which the Earth, the star, and the planet are along 
a straight line. In this configuration, the light reflected by a ring 
system with inner radius r,- and outer radius Tq is, in a simple ap¬ 
proximation, given by the light reflected from a uniform inclined 
disk with radius rg subtracted by the light reflected by a similar 
uniform (and inclined) disk with a radius r,. This reflected sig¬ 
nal also depends on the geometric albedo of the disk/ring system 
(Ap and on the tilt angle (p between the Earth-star-planet line and 
the plane of the rings at the moment of opposition. (The lower 
the value of p, the lower will be the “cross section” of the ring 
as seen by the star and by the observer. By definition, we have 
g{a) = 1 at the maximum phase angle. As such, the total planet- 
to-star flux ratio of a planet with a (optically thick) ring system 
can be approximated by 


^ planet+ring 


lRp\^ 

/ Kn\^ / f'i\^ 

— +A''gr((f>) 

- - - 

\a j ® 

\ a/ \a J 


( 2 ) 


where gr((p) ~ sin^0 is a “reflectivity” function that depends 
on the tilt of the rings with respect to the line of sight (see sec- 
tion [3.1[ l. This model is very simplistic and only serves to under¬ 
stand whether a ring system can explain the observed signal. 

To check that this configuration can explain the observed re¬ 
flected light CCE of 51 Pegb, the first thing we need is to con¬ 
strain the possible values for r, and (i.e., of the inner and outer 
radii of a possible ring system around 51 Pegb). Looking at the 
case of Saturn, we see that r, can be very close to the radius of 
the planet. We see no reason for this to be different in the case of 
51 Pegb, so we thus consider that r, = Rp. 

Prom a dynamical stability point of view, to constrain the 
outer radius we use the Hill approximation. The Hill radius is 
derived from 


77 = a 


Mp 


(3) 


where Mp and M, are the planet and stellar masses, respectively. 
Assuming that the real mass of 51 Pegb is 0.46 Mj^p, that the 


mass of its host, 51 Peg, is 1.04 Mq _ 

that the orbital separation is 0.052 au ( Martins et al.|2015] l, 
then derive a value of H=5.9 R j^p. We s 


Santos et al. 2013[l, and 


we 


lould note, however, that 


several dynamical studies have pointed out that the outer edges 
of the Hill sphere are unstable (see discussion in [Schlichting &[ 
Chang|201 1[ Kenworthy & Mamajek[[20T5|). If we assume that 


only regions within 2/3 H are stable, then the outer edge of the 
ring system around 51 Pegb should be ~4Rjup. 

We note, however, that for radii higher than the Roche radius, 
we should expect that ring particles gather to form satellites. The 
Roche radius, below which a given satellite of density p will 
break up, can be derived from (e.g., de Pater & Lissauer|2010[l 


Rroche — 2.44 Rp 



1/3 


(4) 


where pp is the density of the planet. Assuming that 51 Pegb 
has a radius of 1.2 Ry„p, we derive pp - 0.6 g cmT^ . Considering 
p - 3g cmT^ (typical of rocks), this implies a Roche Radius 
of ~1.5 Ry„p. If we assume that beyond 2 Rroche there should no 
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longer be any ring particle^ this would imply ro„, ~ 3.0Rjup, 
a value lower than the one round above when assuming the Hill 
radius. 

In addition, these values for H and Rroche also depend on the 
real mass for the planet. The values above were computed as¬ 
suming a mass of 0.46 My„p for 51 Pegb. This cotTesponds to an 
orbital inclination of 80 degrees. However, the inclination found 
in Martins et al. is affected by large error bars (i = 80^}°). For 
instance, for values of i = 61 degrees (the lower bound), the real 
mass of 51 Pegb would be 0.53 Mjup. This corresponds to a vari¬ 
ation on the order of 10% in mass, a value that produces a minor 
effect in the derived H and Rroche- 

Adopting values for ro -3 and r, - \ Rjup, as derived above, 
we then estimate a rough value for the expected Fpiamt+ring!F» 
from Eq.|^ (see Fig.[^. To do this, we also need to assume a 
value for Ag and Ag, as well as an inclination (p of the ring system 
with respect to the Earth-star-planet line. We should note that (p 
is at most equal to the inclination of the ring relative to the orbit 
and that this upper limit is only reached at equinox. Eurthermore, 
this problem is highly degenerate. Different combinations of the 
albedos and inclinations will be able to replicate the observed 
(planet -H rings)-to-star flux ratio. 

To keep the different parameters within physically realistic 
values, we decided to fix to 0.3, a value that has been observed 
is several hot-Jupiters (e.g., Cowan & Agol|20TT] [Demory et al. 


2013 I. As an example, using this value for Ag, (p - 60 degrees 
and a value of A^=0.7, we can explain the observed flux ratio 
(assuming rg and r, of 3 and 1 Rjup, respectively). 

Increasing the inclination (p (i.e., increasing the angle be¬ 
tween the rings and the Earth-star-planet line) would imply that 
the projected area of rings would also increase, and lower val¬ 
ues of Ag would be necessary to explain the signal. In Eig.P 


igfi] 

tisfy 


we show the values of the inclination (p against Ag that satisty 
the observed flux ratio. In the hgure, the errors on the detected 
signal-to-star flux ratio were computed from error propagation 
from the recovered values of the signal’s amplitude and EWHM, 

i-e-. f - + where E , Amp, and EWHM 


are the signal-to-star flux ratio, the amplitude, and EWHM of the 
detected signal, respectively. The etTors in the stellar parameters 
were ignored because they are much smaller than the ones of the 
detected signal. Within the three-sigma error bars, possible solu¬ 
tions include pairs of (p and A'g values as low as 0 = 40 degrees 

and of the A^ = 0.6. 
y 

It is not simple to understand what could be reasonable val¬ 
ues for Ag. Observations of Saturn’s rings are not much help in 
this case, since they are rich in ices. We found no discussion in 
the literature about the expected albedo for silicates and other 
refractory species at the equilibrium temperature of 51 Pegb, 
even if such species (e.g., SiOy) are able to condense at the 
equilibrium temperature of 51 Pegb (~1200K - see Pig. 1 in 
[Schlichting & Chang|201 1| ). However, |Draine| ( |l985| l computed 
values for the single scattering albedo of silicates in the inter¬ 
stellar medium that are as high as 0.8 at optical wavelengths. 
These values could be increased if signihcant backscattering oc¬ 
curs near opposition, as seen on Saturn’s rings and other solar 
system bodies (Hameen-Anttila & Pyykko 1972[ Buratti et al. 
|1996[[Wrbiscer et al.|2005|[ 

It is interesting to derive the expected ring rotational velocity 
and compare it to the value of the observed EWHM. Assuming 


^ The rings of Saturn extend beyond the Roche limit, in particular the 
E ring, though they are essentially composed of micron and submicron 
particles (Hedman et al.|201^. 


Keplerian rotatiorjq the velocity of the rings should be 30 km s ' 
at r,- and 17kms“^at Kq. This shows that the presence of rings 
should also signihcantly broaden the CCP of the planet, as ob¬ 
served. We should note, however, that the expected prohle of the 
observed CCP may actually present two different components: 
one produced by the light reflected on the planet disk (that will 
have a EWHM similar to the stellar CCP if the planet rotates 
slowly) and the second produced by light reflected by the ring 
system, very likely with a broader CCP. 


Schlichting & Chang (20111 also call attention to the pos¬ 


sibility that the Poynting-Robertson drag slowly removes dust 
from a ring. As discussed in Burns (19841 and Sfair et aL]! 2009 


and references therein), this effect implies that dust particles wif 
remove orbital angular momentum and spiral into the planet. 
Por a typical hot-Jupiter and assuming a ring orientation of 45 
degree s with respect to the orbital plane, [Schlichting & Chang 
(2011 1 computed ring lifetimes in the range of 10'-10“ years 
(see their Pigure 4). Assuming a lower inclination and a higher 
density optically thick ring, however, this lifetime will increase 
strongly. We thus see no strong reason for a ring system around 
51 Pegb not to resist Poynting-Robertson drag. It is also worth 
noting that, as happens with Saturn, the presence of putative 
shepherd moons in the rings might greatly increase the lifetime 
of the rings; the presence of moons around hot-Jupiters is, how¬ 
ever, questionable ( Weidner & Horne|2010| l. 

These numbers show that the observations reported in 
Martins et al. (2015|l can in principle be explained if we assume 


that 51 Pegb has a ring system. We should note, however, that 
these estimates should be seen mostly as qualitative. Our ma¬ 
jor goal at this stage was to understand if the order of magni¬ 
tude of the effect could be explained using the ring model. Por 
instance, it would be relevant to understand if such a scenario 
should have been detected by existing phase curve observations 
using IR bands ( Cowan et al.||2007) l. The large uncertainties in 
the observed CCP parameters, the properties of a ring system 
(albedos, inclinations), and in the properties of the atmosphere 


^ Vrot - with r the orbital radius of the ring particles and G 

the gravitational constant, valid when the mass of the rings is much less 
than the mass of 51 Peg b. 





Figure 1. Possible conhgurations of the inclination cp and A’g 
given the detected signal-to-star flux ratio. The color gradient 
denotes the 1-, 2-, and 3-sigma error bars. 
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of 51 Pegb (e.g., the albedo and wind velocities) precludes any 
deeper insight into this issue. 


3. Tilted rings? 

In the previous section we assumed that a putative ring system 
around 51 Peg b could have any tilt angle, following the sugges¬ 
tion of Schlichting & Chang ( 2011| l. As can be seen in Fig.[T] 
this has a strong impact on our results. We therefore decided to 
verify this assumption. 

The initial spin state of the planet is unknown. The rotation 
period is supposed to be short, but the obliquity (the angle be¬ 
tween the equator and the orbital plane, here denoted by e) can 
assume any value, due to large impacts and planet-planet scat¬ 
tering at early stages in the formation process (e.g., |Dones &] 
Tremaine]|1993|l. However, due to the proximity of the star, the 


spin of hot-Jupiters slowly evolves until an equilibrium conhg 
uration is reached, corresponding to synchronous rotation and 
zero obliquity (e.g.. Hut 1980|l. The typical time scale 


for 


reaching this hnal equilibrium is given by Correia (e.g., 2009)1 


Port Q 
9nq k 2 


with 

^ Mn\ a 


(5) 


where Port is the orbital period, Q the dissipation quality factor, 
and ^2 the second Love number for the pote ntial. For Jupiter, as - 
trometric observations provide Qlk 2 ~ 10^ ( Lainey et al. 2009 1 . 
Adopting this same value for 51 Peg b gives t ~ 10^ yr, strongly 
suggesting that it has reached its final configuration for the spin. 
The same is true for all hot-Jupiters with a < 0.1 au. Values of O 
as high as 10^ have been proposed for stars (Penev et al. 2012 1 , 
and this value could lead to t ~ 10® yr if we assume ka ~ 0.1, as 
expected for giant planets ( Yoder|1995| l. 

Ring systems are believed to have several possible origins: 
the result of impact events (e.g., Tiscareno|2013 1, captured ob¬ 
jects or satellites that are tidally destroyed (e.g., Charnoz et al.| 
|2009b) Canup||2010 1 , or even remnants from planet formation 


(though this last hypothesis is less likely - Charnoz et al. 2009a i. 


In all cases, they settle in a special plane around the planet, called 
the Laplacian plane (e.g., Lehebel & Tiscareno|2015| l. It is usu¬ 
ally defined as the plane normal to the axis about which the pole 
of a satellite’s orbit precesses ( Laplace||1805| l. For circular or¬ 
bits, the dynamics of the rings’ particles is essentially governed 


by a single parameter, often called the Laplace radius (Tremaine 
let al.|20091l 

Rl = Rp . ( 6 ) 

For To < Rp, the rings can settle in the equatorial plane of the 
planet or in polar orbits. For r, > Rp, the rings can only settle in 
the orbital plane of the system (implying (p - 0). For r, < Rp < 
ro, we expect a transition between the different regimes, called 
“warped” ring. 

The parameter J 2 is related to the flattening of the planet. 
For tidally evolved synchronous planets, we have (e.g., |Correia 
& Rodrlguez|20131l 

J2 = 5k2ql6. (7) 


For Jupiter-like planets k 2 ~ 1/2 (e.g., |Yoder| 19 ^ 5 ) 1 . Using this 
value to compute the J 2 gives Rp - 0.84 Rp for the Laplace ra¬ 
dius (Eq.|^. We thus conclude that for any hot-Jupiter r, > Rp, 
so ring systems can only be observed in the orbital plane (</> = 0). 

This result shows that “warped” rings are unlikely for close- 
in planets such as 51 Peg b, except if we assume that the sys¬ 
tem is young and not yet synchronous. For a coplanar ring sys¬ 
tem, however, according to the approximation presented in Eq.|^ 


we expect no reflected light at all. The approximation is thus 
no longer useful, though it hints that a ring conhguration is 
likely not able to explain the reflected light signal as observed in 
|Martins et al.| ( |2015 1. It is, however, worth understanding what 
is the real amount of reflected light from a ring system in such a 
situation. 


3.1. Reflectivity 

In this section, we no longer assume that the Earth is aligned 
with the star-planet radius vector. It is then necessary to distin¬ 
guish two inclination angles of the rings’ plane: the first one, 
denoted 0,, is computed with respect to the direction of the star, 
while the second, pe, is given relative to the line of sight. We 
note that at conjunction, (pe - (pi - p. Eor a distant star (point 
source), the flux received by the rings depends only on the 
angle pi between the direction of the star and the rings’ plane. 
We thus have 

Fripi) = Fp sin pi , (8) 

where Fp is the flux that the ring would receive if it was per¬ 
pendicular to the incident light. The amount of flux reflected by 
the ring in the direction of the observer also depends on the an¬ 
gle pe between the line of sight and the plane of the ring. If the 
scattering is isotropic, the reflectivity reads (see Appendix) 


qr{.pit pe) 


FriPi) 

Fp 


sin pe - sin pi sin pe . 


(9) 


In the approximation pi - pe - p, as in the previous sec¬ 
tion, we recover the dependency in sin^ p. If a ring system 
around 51 Pegb needs to be coplanar with the orbital plane, 
then sin^ p - 0, and the reflected light flux is thus only due to 
the light reflected by the planet (Eq.j^and Eig.|^. However, for 
hot-Jupiters, the star cannot be seen as a point source, since the 
planet is close enough to receive light coming from the fraction 
of the stellar disk that illuminates the rings. It is thus interesting 
to estimate the total illumination, i.e., a more general expression 
for Frip), and try to understand if this could actually be respon¬ 
sible for the observed signal. 

We let Is be the intensity emitted by the stellar surface and 
we assume that it is uniform, i.e., we neglect the limb-darkening. 
Then, the energy received by a ring element of surface dA is 
given by 


dEripi) = Is 



(ko ■ It) 


i-n ■ k)dA 2 

IIABIP ° 


sin 00 dOo dipo , 


( 10 ) 


where the double integral is computed over the portion of the 
stellar surface visible from a ring element. In this expression, 
ko is the normal of the surface of the star at a given point A of 
spherical coordinates {Rq, (fo, 6 q), B is a point of the ring, and 
n the unit vector along AB (see Eig.|^. Eor 51 Pegb we have 
a - 11 Rq, so we assume that a » Rq. Moreover, as an order of 
magnitude, we only compute the incoming energy at the center 
B of the ring. The flux received by an element of the ring is then 
defined as 

FAPi) = dEAPi)ldA, (11) 


and Fp is equal to nlsfRola)'^. To compute the integral (10 1 , we 
consider two cases. If the ring’s inclination pi is greater man the 
angular radius of the star pc - atan(Ro/a), each element of the 
ring gets the light from the full stellar disk: 


0 < tfio ^ 271, and 0 < 60 ^ tr/2 . (12) 
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Figure 2. Geometry of a system containing a hot-Jupiter with rings, is the unit vector along the radius vector of the planet relative 
to the star, is perpendicular to in the plane dehned by and k {u^ is not necessarily the normal of the orbit). Uy complete the 
orthonormal base frame {u^, Uy, u^). 



Figure 3. Flux received by the ring as a function of the tilt angle. 


In that case, at third order in Ro/a, the reflectivity of the ring is 
still given by gri(f>i, <pe) - sin (pi sin (pe (see Appendix). 

On the other hand, if the inclination 0, of the ring is less 
than (pc, a part of the stellar disk is occulted. In that case, the 
reflectivity becomes (see Appendix) 


^ , sin0ef/7r \ . 2Rq , 


2Rq 


+ cos ipa sin 4>i 


sin ¥> 0-2 'fa Vi 


Ro( . 8 

+ - TT sin -cos (fa sin pi I cos iPa COS pi 

a \ 3 


(13) 


where pa is dehned as 


Rq sin Pa — a tan pi . 


In particular, for small tilt angle {(pi RqIo), we get 


gA^Pi, 0e) 


\37ra 



sini^e ■ 


(14) 


(15) 


For 51 Peg, in the limit of small tilt angles, gr ~ 0.02 sin 0;,, 
that is, the rings only receive about 2% of the maximal Hux com¬ 
puted at pi - 90°. This value is far too small to explain the ob¬ 
served signal as derived in Sect.j^ 


4. Rings from transit surveys 


Even though the results presented above do not support the ring 
hypothesis to explain the signal observed in 51 Peg b, the dynam¬ 
ical discussions presented also show that coplanar rings could be 
present around hot-Jupiters. One can thus wonder if rings are ac¬ 
tually a frequent phenomenon around these sort of planets. 

In a recent paper, [Kenworthy & Mamajek| ( |2015[ l found ev¬ 
idence that the young pre-main-sequence star J1407 (ISWASP 
J 140747.93-394542.6 J1407) may have a planet with a massive 
ring system. This case is not fully comparable to 51 Pegb, in the 
sense that our target is much older and has a much shorter orbital 
period. However, this example shows that present-day photome¬ 
try is able to detect the presence of rings around exoplanets. 

This issue has also also been discussed from a modeling 
point of view (Arnold & Schneider 2004 Barnes & Fortney 


2004 


Dyudina et al.||2005[ Ohta et al. 2009 Tusnski & Valio 


20141. In particular, simulations have shown that if massive 


rings are present in hot-Jupiters, the precision of transit surveys 
like Kepler would have already allowed them to be detected. 
The question is then to understand if we have actually already 
detected rings photometrically but their existence passed unno¬ 
ticed. 

In a recent paper, Zuluaga et al. ( 2015| l have shown that the 
presence of rings would produce (at least) two different effects. 
One of these is that ringed planets would imply that the value 
of the stellar surface gravity (or stellar density) derived from 
the light curve would be systematically smaller than the one 
observed using asteroseismology, for example. Indeed, and al¬ 
though other parameters may be responsible for the derivation 
of erroneous values of the stellar density from the transit light 
curve ( |Kipping|2014| l, such a trend is expected if rings (even if 
not very massive/wide) are common around short-period exo¬ 
planets. 

The other effect discussed in [Zuluaga et al. ( 2015| l is more 
“obvious”: a ringed planet would trigger changes in the transit 
light curve (with respect to a simple planet). In particular, transits 
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Figure 5. Synthetic light curves simulating the light curves of transiting planets with (right) and without rings (left). The parameters 
of each planet are denoted in the inset. Dots correspond to data binned with a 30-minute cadence similar to long cadence of Kepler. 
Different system inclinations are used. In all the simulations, the ring is considered to be coplanar with the orbital plane. 




Impact parameter (b) 



0.0 0.2 0.4 0.6 0.8 1.0 
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Figure 4. Transit depth, duration signal, and derived stellar 
mean density for a transiting planet-i-ring system as a function 
of the impact parameter as computed using SOAP-T-hR and 
EXORING. 


To test these scenarios we modified the SOAP-T tool 
( Oshagh et al.|2013' i in order to add a planetary ring to the tran¬ 
siting planet (hereafter we call this code SOAP-T-hR). SOAP-T 
was originally designed to generate the radial velocity variations 
and light curves for systems consisting of a rotating spotted star 
with a transiting planet. The model assumes that the rings are 
uniform and completely opaque and that they have an orienta¬ 
tion with respect to the orbital plane. By comparing the transit 
light curves of SOAP-T-hR with those of SOAP-T, we are able 
to recognize the impact of rings on the transit light curves. A full 
description of this code is beyond the scope of the present paper. 

To check that the SOAP-T-i-R code was working properly, 
we compared its results with the ones obtained with the avail¬ 
able EXORING code ( Zuluaga et al.|2015 i. The obtained transit 
duration and depth are shown in Fig.as simulated using both 
codes as functions of the impact parameter for a ring system 
that is coplanar with the orbital plane. The comparison of the 
results obtained using SOAP-T-i-R and EXORING for the tran¬ 
sit duration show very good agreement. On the other hand, the 
transit depth obtained from SOAP-T-i-R displays deeper transits 
than those using EXORING. This difference is explained by the 
fact that the stellar limb darkening is neglected in the EXORING 
code, while in SOAP-T-i-R we consider a quadratic limb darken¬ 
ing coefficient close to the solar value. Indeed, if we assume no 
limb darkening in SOAP-T-hR, we obtain the same results as 
EXORING (see also Fig.|gi. 

In Fig. 1^ we show the results of our simulations after com¬ 
paring the transit light curves produced by a planet (left) and by 
a planet with rings (right). Here we assume a ring system that 
is coplanar with the orbital plane and that has inner and outer 
radii of 1 and 4Ryj,p, respectively. Quadratic limb darkening pa¬ 
rameters ul=0.29 and u2=0.34 were used in this simulation (as 
expected for a Sun-like star). As seen in the figure, the impact 
of such a ring system can be quite significant. The most inter¬ 
esting result is that the shape of transit light curves of ringed 
planets (deep transit, long duration, and shape - in most cases 
“V” shaped) look very much like the eclipse light curve of one 
eclipsing binary. Therefore, possible transiting planets with rings 
could have been observed by transit surveys, such as Kepler and 
CoRoT, however they could have easily been misclassified as 
false positive candidates. 


should be deeper and longer. If the large rings are present, the 
amplitude of the transiting signal could even be similar to the 
one expected for an eclipsing binary star. 


We need to add, however, that the actual capability to distin¬ 
guish between an eclipsing binary and a transiting ringed planet 
needs to be assessed by, for example, simulating the expected 
light curves in detail (including the different sources of noise) 
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and investigating the residuals when fitting both models to the 
simulated data. We leave this detailed analysis to future studies. 

The lower panel of Fig. [^presents the impact of the “unac¬ 
counted” effect of a ringed planet on the derived stellar density 
as a function of planet impact parameter. To estimate this we 


used the stellar density as derived using the Eq. 9 in Seager & 
Mallen-Ornelas (2003|l. The plot shows that the stellar density is 


underestimated as we move toward higher impact parameter val¬ 
ues. In this respect it is interesting to note that Huber et ari ( |2013| ) 
find that the difference between the light-curve stellar density 
and the value derived using asteroseismology is a function of 
the impact parameter of the planet. We are not advocating, how¬ 
ever, that rings are the definite explanation for this trend. In any 
case, most of the systems in the Huber et al. paper are low-mass, 
small-radius, and not Jovian in nature. 


5. Conclusions 


In this paper we have explored the possibility that the reflected 
light spectrum observations of 51 Peg b ( [Martins et al. 2015| l can 
be explained if we assume that this hot-Jupiter has a ring system. 
Using a simple model we showed that, overall, the observed sig¬ 
nal can indeed be explained by a ring system under the assump¬ 
tion that rings are tilted with respect to the orbital plane of the 
planet. We showed, however, that dynamical arguments suggest 
that in any synchronous hot-Jupiter like (we expect) 51 Pegb, 
this configuration is unlikely. In the case of a ring system copla- 
nar with the orbital plane of the planet, we also showed that the 
total amount of incident flux is about two orders of magnitude 
smaller than the one needed to explain the observations. 

The study shows, however, that the analysis of the reflected 
light spectrum from an exoplanet could be a very interesting 
method of detecting rings around short-period systems, in partic¬ 
ular those not transiting. This approach can also complement the 
measurement of brightness variations along a phase curve, as al¬ 
ready proposed by Dyudina et al. ( 2005| l. Observations of appar¬ 
ently high-albedo planets that present broadened spectral lines 
or line profiles showing two different components (see Sect.|^ 
could hint at the presence of rings around other planets. A de¬ 
tailed model of the observed line-shapes could indeed provide 
relevant information about the system. 

We also discussed the possibility that planets with rings 
could have been detected by space missions like Kepler but 
simply discarded as binaries owing to the shape and depth of 
the transiting signal. In this respect we propose that it would 
be interesting to obtain precise radial velocity measurements of 
candidate binary stars from the Kepler field, in order to derive 
the masses of the companions. The study of the light curves to 
search for transiting secondary binary like-signals with no sig¬ 
nature of secondary eclipses or beaming and ellipsoidal effects 
( |Mazeh et al.|2012) could also help to select the best candidates. 
We note that smaller ring systems could also be responsible for 
slightly deeper transits, leading to deriving an inflated radius for 
the transiting planet. A careful analysis of the data could be rel¬ 
evant in such cases. 
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For this problem, we set 


a tan <^>1 




where 
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The integrand in the expression of Fr is of order unity. We ex¬ 
pand it at the first order in Rq! a ^ 1. We get 


Fr^L 


Rl 


JJ' ^ sin (pi cos 00 H—- [ sin (pi{3 cos^ 0o - 1) 

-H COS (pi sin 00 cos 0o cos ^o] j sin 0o £/0o ■ (A.3) 

We consider the most general case where each element of the 
ring only sees a fraction of the stellar disk (see Figure [ATT) . This 
case happens when the tilt angle <pi is less than the angular radius 
of the star (pc - atan(/?o/a)- In this configuration, the visible 
surface is delimited by two curves: the arc CDE in the 
of the star and bounded by -njl - ipa ^I + Pa, and the 

arc EC, which is half of a circle of radius cos ipa in the plane 
of the ring. For commodity, we recall the definition of the angle 
Pa given in Eq. ( [T4| 

Rq sin Pa - a tan (pi . 


ff VxA-dL^S A-dt. 

JJl. JdS 
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where 
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1 Rq 
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R^O. 

For the line integral CDE, we use r = (x, y, z) with 
X = cos tpo , 0 = sin i/jQ , z = 0, 


(A.6) 


(A.7) 


Figure A.l. Area enclosed by the thick curve CDEC is the visi¬ 
ble part of the star seen from an element of the ring. 
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Appendix A: Flux received by the ring 

Here, we detail the computation of the flux Fr{(pi) received by 
the ring per unit area. The notation is the same as in Sect. 3.1 of 
the main text (see also Eigurej^. The general expression of the 
flux derived from Eqs. ([T§and([^ is 

32 


where (fQ goes from {-njl — pa) to (;7r/2 -i- (fa),- While for the line 
integral EC, we set r - (x, y, z) with 


X = cos (fa sin (pi sin xp — sin ipa , 
y — cos ipa cos (p , 

Z — cos iPa COS (pi sin Ip , 

where ^ ranges from 0 to n. As a result, we get 


(A.8) 
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In the case where 0, > (pc, Eq. ( |A.9| l still holds if we set pa - nPl 
so we get 

IR 

F,(0,) = 7r/,^-^j sin0,', 
while for (pi « (pc, Pa ~ a(pilRo and 

Fc{(Pi)^Is[a] 3 ^ + 2 <^' 


Appendix B: Reflectivity 


The reflectivity of the ring is computed by assuming an isotropic 
scattering. Eurthermore, it is assumed that given an incoming 
flux Fy, only a fraction is re-emitted in the visible spec¬ 
trum. Thus, the luminous intensity C of the rings is uniform and 
such that 

A[Fr((Pi)^7lIc. (B.l) 

Besides this, the flux received on Earth from the disk is 


Fringitpi, (pe} — 


’ proj 


(0c) 


D2 


(B.2) 


To compute the surface integral ( |A.3| l, we make use of the 
Stockes theorem that transforms a surface integral over E into a 
closed integral over its boundary dl, as 


where the ratio of the projected surface 5pioj of the rings on the 
plane of the sky divided by the square of the distance D to the 
Earth represents the solid angle under which the rings are seen. 
As a result, we get 


(A.4) 
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Moreover, the stellar flux F* received on Earth is 




a Y 




(B.4) 


where we used Fp - nlsiRolar. Combining Eq. (B.3 1 and (B.4 1 , 
we get 


“ 0e) 
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with gr a function representing the reflectivity of the rings given 
by 

FJ(pi) 

gricpi, (t>e) = ^ Sin ■ (B.6) 
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